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THE ANALYST. 

Vol. VI. July, 1879. No. 4. 

REDUCTION OF OBSERVATION EQUATIONS WHICH CON- 
TAIN MORE THAN ONE OBSERVED QUANTITY. 



BY CHAS. H. KUMMELL, ASST. U. S. LAKE SURVEY, DETROIT, MICH. 

The ordinary case, in the application of the method of Least Squares, 
where a quantity is observed which is a known function of certain unknown 
quantities, and the case, occurring principally in higher geodesy, where the 
observed quantities have to fulfill rigorous conditions, have received by 
Gauss, Bessel, Hansen, Encke and others such an exhaustive treatment, 
that it would be difficult to add anything of importance to their labors. 
The case mentioned in the heading of this article, which is nothing but an 
extension of the ordinary case, whenever it occurs, is generally forced into 
the treatment of the ordinary case, either by considering one class of observ- 
ed quantities without error (comparing standards of length at different 
temperatures, we really have two observed quantities, viz., difference of 
length and temperature, but in this case the temperature is considered in the 
reduction without error), or some function of the observed quantities is treat- 
ed as an observed quantity. — In Mr. Hill's reduction of Eegnaults' experi- 
ments on the volume of atmospheric air, Analyst, Vol. IV, No. 4, there 
are three classes of observed quantities, viz., Vol., V, Pressure, P and tem- 
perature, T. In the first class of experiments, where the temperature is 
nearly the same there are then only two variable observed quantities, V 
and P, and in applying Least Squares Mr. Hill makes the sum of the 
squares of the corrections to log(PF) a minimum. — Upon close examina- 
tion it will be found that a great many observations are of a more or less 
complex nature so that the observation equations contain more than one ob- 
served quantity, and it may be of interest to show how such observation 
equations should be treated. 
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Suppose we have the following m observation equations : 
f x {a, b,... »!+&>!, y 1 +Sy l , . . .) = 0, 
/ 2 (o, b, . . . a: 2 + <Ja; 2 , y 2 + dy 2 , . . .) = 0, 



(li) 
(1 2 ) 



/ m (o, b, . . . a? ra +^ar m , y m +dy m , . . .) = 0, (1„) 

where the form of these functions may be considered different, provided the 

same unknown quantities, a, b, either all or in part, occur in them. 

Also the observed quantities need not be of the same class in all observation 
equations but I suppose here each observed value to enter only one of them. 
Since the corrections dx lt Sy 1} . . . §x 2) dy 2 , . . . must be small quantities, 
we have, iff lf f 2 , . . ■/„ denote the values of (l^, (1 2 ), .. . (l m ), if the un- 
corrected values, x t , y x , . . . x 2 , y 2 , . . . x m , y m , ... are substituted, 

/-+g*^*.+ ■■■ = «• ay 

Let p x = weight of » 1 ; g^ = weight of y t ; . . . (2 t ) 

P2 = " " ** ; ?» = " " y. ; • • • (2 2 ) 

J>„= " " * m ; 9m - « « y.;... (2.) 
The quantities/^, / 2 , . . . .f m are independent with respect to the obser- 
vations as we here suppose. Let P lf P 2 , P 3 , P m be their weights 

respectively then we have, evidently, 

p.-pMxJ + qi \dfJ + '" (3l) 

p~ 2 = p~ 2 \^) +tM?J + * " (3a) 

which expressions are in general functions of o, b, . . . and of the observed 
quantities. 

According to the principle of least squares we must make 

[pdx 1 ] +[?%*]+ . . . = a minimum, (4) 
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or, by (I') [pfe>]+ [g( /+(<?/ ^ ) ^ + •• •)'] + • • • = a min. (4') 

Equating to zero the differential coefficients with respect to dx ly dx 2) . . . 
dx m , successively, solving for these corrections, and doing the same for Sy lf 
8y 2 , . . . dy m , etc., we obtain regarding (3) 

& = _AA . ^A . dy = _i±k ,iL (5 ) 

1 Pi <V yi q x d Vl 

8x 2 = -iik. . fk; % 2 = -lik . d Ji (5 2 ) 

Pz efe 2 J " q 2 dy 2 



X„ " m fm df m # jj "mfm df m 

Pm dx m J q m dy m 


(5 m ) 


and substituting these expressions into (4) we obtain 




[P/ 2 ] = a minimum ; 


(4") 


4£v\\=°> 


(«-) 


iM=°- 


(6.) 



Only in the simplest case these conditions may be employed for the direct 
solution of our problem. Before attempting the general problem, I take 
such a simple case. Let 

f 1 = ax 1 —y 1 +b = —adx 1 +dy 1} (a x ) 

f 2 = ax 2 —y 2 +b = —aSx 2 +dy 2 , (o 2 ) 



f m = ax m — y m +b = —adx m +8y m , (a m ) 

be the observation equations, then we have by (3) 

1 a 2 . 1 th \ 

r 2 Pa ?2 



-p- = —+—- ( 6 «) 



by (4") ~^r m = 

L y q J 



"We have then by (4") a 2 , 1 = a minimum,* 



*If p = q = 1 we have [(ax — y+6) 2 -n(a 2 +l)] = a minimum. This is Mr. Adcock's 
principle of making the normals as short as possible. See Analyst, Vol. V, p. 54. 



-100- 



and by (6) 



d 
2da 

A 
~2db 



- (ax—y+b)* 

L p g 

(a x — y-{-b) 2 
~~ a* 1 

P 2 



■]-. 



or 



[fd "*-*+" )] - [W&( —*+»)'] - 0> 



+ 1 ; 



aa; 



-y + 6 






+ 1 
9 



]... 



(C) 



(«») 



K) 



W 



Unless we make a certain assumption with regard to tlie weights of the 
observed quantities, we cannot solve, by a direct process, even these eq'ns. 
Such an assumption is : that the precisions of x and y which enter the same 
observation equation, are considered to have a constant ratio. The x and y 
observations having been made as nearly as possible simultaneously, they 
are made under the same conditions and are therefore as good as those con- 
ditions permit. From this however it does not follow that they have equal 
weight, more particularly if they are heterogeneous quantities, as for exam- 
ple, temperature readings and microscope readings on the comparator in 
comparing standards of length, or signalling the transit of a star over a 
thread in a telescope on the chronograph and reading the level, in time ob- 
servations, etc. If weight is defined as the reciprocal square of the probable 
error, it is clear that the unit in which a quantity is expressed makes a great 
difference in the weight to be assigned. All we can safely assume therefore 
with regard to the weights of simultaneous observations that enter one 
observation equation is this : that there is, throughout, a constant ratio in 
their precisions, the square of which is the constant ratio of their weights. 
To obtain this we may take the reciprocal ratio of the practical limits of 
error in each class of observations for the ratio of the precisions. Such 
limits are known sufficiently near by any experienced observer. Even a 
computer can draw some inference from the results given to him, which 
will enable him to get some rough idea of this ratio ; for a good observer 
will not give his results with any more accuracy than he thinks he can es- 
timate. Now suppose x to be given to p., and y to v decimal places, what- 
ever be the units in which they are expressed; then it is the computer's 
duty not to alter the respective figures (not the absolute values) any more 
than necessary to establish consistency. This he will effect if he takes y 
10 v — .« times as precise as x. If however, for example, y is the mean of It 
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directly observed values, he should take y 10"— ^h times as precise as x, 
and similarly he might reason in other more complex cases. Suppose then 

9i = k Pi, («i) 

?» = *Pa> («a) 

q m = ty m , (e m ) 

then we have, substituting in (d„) and (<7 6 ), 

a[px !r \—[pxy]+b[px']— - r _g— J [p/]— 2a[pa^]+a*[pa; 2 ]— 26[>y] 

+2a6[ pa! ]+6 J [p] }=0, (««) 

a[jw]-[py]+6[p]=0. (<?) 

From (dj) we have & =(l^[p])*J [py] — a[px] )■ (f) 

Substituting this into (d^), we obtain, after reduction, 

+a \~ ; [po;y]-(l-f-[p])M[p2/] f £~ U ' W 

and solving we obtain 

"- 2Qwy.l] -\U + ~~ 4[p^.l] a """ J' W 

where use has been made of Gauss' notation of auxiliaries as used in his 
method of substitution. 

To decide which of the signs must be taken for the minimum minimorum, 
I take two extreme cases: 

1. Let the weight of the x be infinite, or let us suppose the x without 
error, then k = in order that q=kp=a, finite quantity. Then (h) becomes 

a' = [ qxy -]} or = oo. (A') 

Iqx 2 . 1J 

2. Let the y be without error, or h = oo, then (h) becomes 

a"= fMl4_] or = 0. Qi") 

[pxy.l] 

Only the results given first, which correspond to the upper sign, can 
be used ; for they agree with the results from the ordinary treatment if 1, 
the (y)'s and 2, the (a:)'s are corrected. "We have then to adopt the upper 
sign in (h). 

Ifp x = p% = . . . .p m = 1 and also k = 1 we have Mr. Adcock's prob- 
lem, Vol. V, page 54, and (h) becomes 

„ - M-(i^-m)|>]'-|y] + Mj' 

2[>y]-(2-*-m)|>][y] 
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+ N i+ 4\ ~ [ay]-(l+»)M[y] ' J l ° ; 

This however is not identical with Mr. Adcock's formula in consequence of 
an error which he committed in eliminating 6. 

Having found a and 6 we compute the / and P by formulae (6) and (c), 
and the minimum sum of squares by (4"). 

We have then the probable error of an observation of weight 1, since 
there are two unknown quantities, 

Although the probable errors of a and b could be deduced in terms of 
the probable error of observation from their expressions just given, yet since 
this requires tedious operations and since the approximate solution which I 
am going to give of the general problem also gives the probable errors of 
the unknown quantities, I pursue this special case no farther. 
Assume 

1*1 1 

•Pi = jy^ • ••(7]); Pi = -jy C^); • • • J 3 - — ~2j^> ( m > 

then (4") becomes [/ 2 -^-^] = a minimum. (4"') 

Let a , b , . . . be some approximate values of the unknown quantities, 
and Ja , Jb , . . . their corrections, then we have, in (4'") exact to quanti- 
ties of the second order, 

K]=K].+I.K]>h-ik]^ + ... 

+ (4™) 

= a minimum ; 
and the conditions for a minimum, or the normal equations, are 

d r/H - d r/H j- d2 r/H a _l d2 r/n ^ 

LIU ~2daLxJo + "2d^L^Jo 0+ 2»LlUo ° 

+ ... = 0, (£) 

K] -i[S.+5»[ j 9. i - + w[ j S].^ 

+ ... = 0, (61) 
*This is done because the N's are easier to differentiate than the P's. 



2Ja 



d 

2J6 
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We have 



2daLiVJo _ LMda 2iT da /Jo' W 

2d6 [_Nj ~ \_~N\db 2N db I J o' W 

rf 2 r/n _ri-^-Z^V^-i-^l 

2dad&|_.#'Jo LN\da N'dalKdb N'dbJJo 

, \~f(d*f_f gga , 8 v 

^LM<fad6 2N'dadb) Jo' y * } 



2db* 



r/n _rite£_/ fflH +yf( d *f-f d2N \~] (g \ 

[_NJo~~LN\db N'dbl J^^NXdb* 2JV' d& 2 >/ J A **' 



Having computed the first and second differential coefficient of/ and N, 
these coefficients can be found and used for the normal equations (6'). 
If we place 

dfi_ /i gjj _ y . <?/i /l^-F /o >, 

da 2N/ da x ' dfc 2N X " di 1 '"" <- lj 

d / 2 _ /j dN 2 _ y • d ^_ /» dN 2 — F CQ ^ 

^a "2^7 da ~~ 2 ' d6 2N 2 , "aZT _ 2 '"- ^ 

djm Jm aJN OT y t aj m j m aJN m y (a \ 

da 2NJ da "' db 2N m 'd6 m '"" l ^ 
we have sufficiently near (replacing now 1-f-N by P) 

[PX/] +[PX 2 ] Ja +[PXF] J6 + ... = 0, (61') 

[PF/] +[PXF] Ja +[PF' ] ^ + . . . = 0, (6'i) 



These normal equations we might consider to have been formed from the 
following fictitious observation equations : 

Jj_ =/i +X 1 Ja + Y l Jb + .... (weight = P x ) (10 x ) 

J 2 =/ 2 +X 2 Ja + F 2 J6 + .... (weight = P 2 ) (10 2 ) 



4» = / M +X> + F »^o + (weight = P„) (10„) 

if [PJ 2 ] is made a minimum. 
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The equations (6"), (6") . . . not only give the corrections Ja , db Q , . . . 
to the assumed values of the unknown quantities, but also their weights in 
the usual manner. 

If n= number of unknown quantities there are m — n conditions between 
the observed quantities alone. We have then for the probable error of an 
observation of weight 1 

where [P/ 2 ] is to be computed by means of the final values of a, b, . . . 
and if p a , p b) . . . are the weights of a, b, ... as deduced from (6"), (6") 
... we have their probable errors, 

r a = -J-, (12.) 

n= " (12.) 

VPt 



Applying this process to the above simple case we have 

da Xl ' db~ ' da* ~ dadb~ db* ' ' " 3l) 

da~ 2 ' ~db~ ' ~da*~~cfadb~~dW~ ' {i2> 



dTS t _2a dN i _ ft d 2 N t _ 2_ d*TS t _ d*N 1 _ . . 

da Pl ' db ' rfa 2 Pl ' dadb " di 2 " ' {%) 

da p 2 '' db ' da 2 p 2 '' dadb db 2 ' {z > 
We have then 



55LNJ, LW N- Pl )jo~L^ / r ~T a ) lo {a) 

* r/ 2 i _ r / i _ r P / 1 / M 

2a^LNJo~L^r "7" / Jo L7 ; Jo' (oo) 

<* 2 r/ 2 i _r P L_2P/ a \ i „ n 

wLU.-L r ~y ) Jo' (o6) 



2d6 2 L Nj L J 



o 
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Substituting into (6') we have 

+ Qp (xJ^f a )~^ Ab, = 0, (m„) 

[p/] o+ [p(,-^/a)] / a 0+ [p] / 6 =0. K ) 
Employing (9) (6") and (10) we have also 

X x =x x -l±3-±a; Y 1 =l, (» t ) 

ri 

x 2 = * 2 -^Aa ; ^2 = 1, (« a ) 

1*2 

Should the results of these two methods differ materially, the indication 
is that there is considerable uncertainty about the determination ot the un- 
known quantities or that their probable errors are large. It would not be 
judicious in that case to attempt a more exact solution. 



SYMMETRICAL FUNCTIONS OF TEE SINES OF THE 
ANGLES INCLUDED IN THE EXPRESSION 

«o+ 

n 



Denoting these angles by a , a lt a 2 , . . . a„_i, the values of sin na and 
cos ni will be the same for all ; hence, putting x for sin d and na for nO in 
the developments of cos nO and sin nO in terms of sin (quoted on page 17) 
we have the eq'ns whose roots are sin a , sin a lf etc. ; viz., when n is even, 

cos na=l _£ *+ ^1 * . . . +(-!)*. ^ 2 )(^. . . KMf_;, 

(1) 
and when n is odd, 

sin n«™ -"*££ * . . . ^^V ^ W K -(^, 

<3! « ! 

( 2 ) 
On dividing by the coefficient of x n each of these eq'ns takes the form 

_ ~» n(n—l) n-2_i_ n(n—l)(n—2)(n—S) ¥ >- A _ 

n 2_( n _2) 2 "•" [ n »_( n _2) 2 ] [ n 2._(n_4) 2 ] ' " " ' 

the equations differing only in the last terms. Reducing the coefficients, 



